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ABSTRACT 

We report a detection of galaxy-QSO cross-correlation wqq in the Sloan Digital Sky Survey (SDSS) 
Early Data Release (EDR) over 0.2 — 30 arc-minute scales. We cross-correlate galaxy samples of different 
mean depths r' = 19 — 22 (~zc = 0.15 — 0.35) with the main QSO population (i'q < 19.1) at ~zq ~ 1.6. 
We find significant positive correlation in all cases except for the faintest QSOs, as expected if the signal 
were due to weak lensing magnification. The amplitude of the signal on arc-minute scales is about 20% 
at ~zq = 0.15 decreasing to 10% at ~zq — 0.35. This is a few times larger than currently expected from 
weak lensing in the ACDM models but confirms, at a higher significance, previous measurements by 
several groups. When compared to the galaxy-galaxy correlation wee, & weak lensing interpretation 
indicates a strong and steep non-linear amplitude for the underlaying matter fluctuations: a ~ 400 on 
scales of 0.2 Mpc/h, in contradiction with non-linear modeling of ACDM fluctuations. We also detect 
a normalized skewness (galaxy-galaxy-QSO correlation) of S3 ~ 21 ± 6 at z ~ 0.15 (S3 ~ 14 ± 4 at 
z ~ 0.35), which several sigma low as compared to standard ACDM expectations. These observational 
trends can be reconciled with lensing in a flat A universe with og — 1, provided the linear spectrum 
is steeper (n ~ 1) than in the ACDM model on small (cluster) scales. Under this interpretation, the 
galaxy distribution traces the matter variance with an amplitude that is 100 times smaller: ie galaxies 
are anti-bias with b ~ 0.1 on small scales, increasing to b ~ 1 at ~ 10 Mpc/h. 

Subject headings: galaxies: clusters: general — large-scale structure of universe — cosmology: 
observations — gravitational lensing 



1. INTRODUCTION 

Weak gravitational lensing by foreground large-scale 
matter density fluctuations could introduce significant 
density variations in flux-limited samples of high redshift 
objects, such as QSOs. This is sometimes called magni- 
fication bias or cosmic magnification. In principle, it is 
possible to separate the intrinsic density fluctuations in 
QSO samples from the weak lensing magnification signal 
by cross-correlating the QSOs with a low redshift galaxy 
sample. This allows a direct measurement of how the 
galaxy distribution traces the underlaying mass distribu- 
tion. Much work have been done in this direction both on 
theory and observations (for recent reviews see Norman 
& Williams 2000, Bartelmann & Schneider 2001, Benitez 
et al. 2001, Guimaraes, van de Bruck and Brandenberger 
2001). 

Correlation of low redshift galaxies and high-redshift 
AGNs or QSOs typically find significant excesses of fore- 
ground objects around the QSO positions (Seldner & Pee- 
bles 1979, Tyson 1986; Fugmann 1988,1990; Hammer & Le 
Fevre 1990; Hintzen et al. 1991; Drinkwater et al. 1992; 
Thomas et al. 1995; Bartelmann & Schneider 1993, 1994; 
Bartsch, Schneider, & Bartelmann 1997; Seitz & Schnei- 
der; Bem'tez et al. 1995; Bem'tez & Martinez-Gonzalez 
1995, 1997; Bem'tez, Martinez-Gonzalez & Martm-Mirones 
1997;Norman & Impey 1999;Norman & Williams 2000). 
See also Burbidge et al. 1990 (and references therein) for 
individual associations and different statistical indications. 
On the other hand, there has also been reports of sig- 
nificant anti-correlation between galaxy groups and faint 
QSOs (see Croom & Shanks 1999 and references therein) 



which could also be explained by weak lensing magnifi- 
cation because of the flater slope of the faint QSO num- 
ber counts. Even though the shape of the galaxy-QSO 
cross-correlation has not yet been well constrained, these 
results seem qualitatively in agreement with the magnifi- 
cation bias effect, but the amplitude of the correlation is 
found to be higher than that expected from gravitational 
lensing models based on ACDM. 

What is the origin of this discrepancy? Part of the prob- 
lem could be due to the lack of well-defined and homoge- 
neous samples. After all we are looking for a small effect, 
may be as low as 1%, and any systematics in the sample 
definition is likely to introduce cross-correlations at this 
level (see §2 below). On the other hand we know little 
about clustering of dark matter on sub-megaparsec scales. 
Could the discrepancy be due to real deviations from the 
ACDM paradigm? 

Recently, Menard & Bartelmann (2002) and Menard, 
Bartelmann & Mellier (2002) have explored the interest 
of the SDSS to cross-correlate foreground galaxies with 
background QSO. The SDSS collaboration made an early 
data release (EDR) publicly available on June 2001. The 
EDR includes around a million galaxies and 4000 QSOs 
distributed within a narrow strip of 2.5 degrees across the 
equator (see Stoughton et al 2002 for details). As the 
strip crosses the galactic plane, the data is divided into 
two separate sets in the North and South Galactic caps. 
The SDSS collaboration has presented a series of analy- 
sis (Zehavi et al. 2002, Blanton et al. 2002, Scranton 
et a.l 2002, Connolly et al. 2002, Dodelson et al. 2002, 
Tegmark et al. 2002, Szalay et al. 2002, Szapudi 2002) of 
large scale angular clustering on the North Galactic strip, 



2 



which contains data with the best seeing conditions in the 
EDR. Gaztahaga (2002a, 2002b) presented a first study 
of bright (</ ~ 20) SDSS galaxies in both the South and 
North Galactic EDR strip, centering the analysis on the 
comparison of clustering to the APM galaxy Survey (Mad- 
dox et al. 1990). 

In this paper we will follow closely Menard, Bartel- 
mann & Mellier (2002, MBM02 from now on) proposal 
to study the galaxy-QSO cross-correlation signal in the 
EDR/SDSS. The paper is organized as follows. In section 
§2 we present the samples used and the galaxy and QSO 
selection. Section §3 is dedicated to the extinction con- 
tamination. Section §4 and §5 presents the main results 
and its interpretation, while §6 and §7 are dedicated to a 
discussion and a listing of conclusions. 

2. QSO AND GALAXY SAMPLES 

The galaxy samples are obtained from the EDR and 
converted into pixel maps of different resolutions as de- 
scribed in Gaztanaga (2002a, 2002b). We select objects 
from an equatorial SGC (South Galactic Cap) strip 2.5 
wide (—1.25 < DEC < 1.25 degrees.) and 66 deg. long 
(351 < RA < 56 deg.), which will be called EDR/S, and 
also from a similar NGC (North Galactic Cap) 2.5 wide 
and 91 deg. long (145 < RA < 236 deg.), which will 
be called EDR/N. These strips (SDSS numbers 82N/82S 
and 10N/10S) correspond to some of the first runs of the 
early commissioning data (runs 94/125 and 752/756) and 
have variable seeing conditions. Runs 752 and 125 are 
the worst with regions where the seeing fluctuates above 
2" . Runs 756 and 94 are better, but still have seeing fluc- 
tuations of a few tenths of arc-second within scales of a 
few degrees. These seeing conditions could introduce large 
scale gradients because of the corresponding variations in 
the photometric reduction (eg star-galaxy separation) that 
could manifest as large scale number density gradients (see 
Scranton et al 2001 for a detail account of these effects). 
We will test our results against the possible effects of see- 
ing variations, by using a seeing mask (see §3.1). 

Redshift targets in the SDSS are selected in r' for galax- 
ies and %' for QSO. Here we use i' for QSO and both i' 
and r' for galaxies. 1 The first choice has the interest that 
both galaxies and QSO come from the same photometric 
reduction and are therefore subject to similar systematics. 
Galactic extinction is also smaller in this band. The sec- 
ond choice provides a comparison with previous results on 
galaxy clustering. To maximize the number of galaxies, 
and therefore the possibility of lensing, we choose broad 
magnitude bands. We will focus our results in comparing a 
nearby sample r' = 19— 17.5 (~z ~ 0.15) with a distant one 
i',r' = 22-17.5 (z ~ 0.35). Galaxies selected with i' < 22 
are almost the same (within few percent) to galaxies in 
r' < 22 as the K-correction cancels out with the color 
evolution (see Fukugita et al. 1996). Here we take the 
i' < 22 band as our nominal choice to minimize extinc- 
tion (this will give larger area coverage, see below). As 
we will show, both galaxy-galaxy and galaxy-QSO cross- 
correlations turn out to be almost identical in i' < 22 and 
r' < 22, the only difference being slightly smaller errors 

1 We will use z' , i', r', <?', u' for 'raw', uncorrected magnitudes, and z 
cording to Schlcgcl ct al. (1998) r' = 18 corresponds roughly to an avc 
E(B -V)~ 0.03. 

2 Recall that i' are 'raw' magnitudes. The extinction corrected limit i 



in the i' < 22 sample. For the bright sample we stick to 
r' < 19 to provide a more direct comparison with previ- 
ous results on galaxy-galaxy clustering. Using the redshift 
distributions in Dodelson et al. (2002), the mean red- 
shift for 17.5 < r> < 19 (with ~ 81000/117000 galaxies) is 
z ~ 0.15, while 17.5 < i',r' < 22 (with ~ 480000 galax- 
ies within our mask) have ~z ~ 0.35. We will center our 
analysis over these two samples, which will be sometimes 
refered to as r' < 19 and i' < 22. 

The 1st and 3rd slices in Figure 1 shows the EDR/S 
and EDR/N pixel maps for i' < 18.5 and 7 arc-minute 
resolution. Note the "barrel" shape in the EDR/N. As 
far as we have been able to find out, this seems an unre- 
ported artifact in the EDR redshift sample release which 
does not seem to include rcdshifts for secondary targets 
in this equatorial strip (see Stoughton et al 2002 for de- 
tails). This is not a problem in our analysis other than we 
are missing a good fraction of the EDR/N because there 
are no matching rcdshifts for the QSOs. The top slice in 
Figure 2 shows a zoom over the central region of EDR/N. 

We recovered the QSO sample form the SDSS/EDR 
as described in detailed in Schneider et al (2002, see 
also Richards et al. 2002). We use point-spread func- 
tion magnitudes for QSO (which will be quoted with 
a subindex: i'g) and Petrosian magnitudes for galaxies 
(which will be quoted without any subindex: ie i'), as 
indicated by above references. We have considered two 
QSO samples. The one obtained directly from the EDR 
data archive (ie just with specobj.specClass = 3) con- 
taining 4275 QSOs, which we call EDR/QSO, and the 
corrected public QSO sample presented in Schneider et 
al (2002), which we called SDSS/QSO and contains 3814 
QSOs. Almost all SDSS/QSO is contained in EDR/QSO. 
SDSS/QSO contains a handful of additional BALQSOs 
and some quasars that were identified by eye during the 
checking process. SDSS/QSO excludes most of the low- 
redshift AGN (Mj < —23) and a small number of narrow- 
line QSO included EDR/QSO, but this makes little differ- 
ent for our sample which has z > 0.8. 

Out of the EDR/QSO sample we impose two further 
cuts: specobj.zStatus > 2 to exclude failed or inconsis- 
tent redshifts and specobj.zConf > 0.7 to exclude red- 
shifts with confidence less than 70% (see Stoughton et al. 
2002 for details). In both samples we restrict our analysis 
to 0.8 < z < 2.5 where the QSO distribution is com- 
pact and homogeneous (the lower cut avoids overlap with 
the galaxy sample). After these cuts the main difference 
between EDR/QSO and SDSS/QSO are the few missing 
radio QSO and a few narrow-band emission QSO's. Over 
97% of the QSOs are the same and have the same param- 
eters. Note that according to Schneider et al. 2002, QSO 
selection criteria were changed systematically during the 
EDR runs resulting in large variation on the number of 
QSOs per field. Nevertheless both samples are photomet- 
rically completed to better than 90% in 16.5 < i'g < 19. 1 2 
(Schneider et al. 2002, Stoughton et al. 2002). We will 
therefore restrict our analysis to 16.5 < i'q < 19.1 and 
further compare results for the two versions of the target 

*, i* , r* , g* , u* for extinction corrected magnitudes. For example, ac- 
rage extinction corrected r* ~ 17.9 for a mean differential extinction 

si* < 19.0. 
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selection that are relevant to our sky coverage: v2.2a and 
v2.7 (see Stoughton et al. 2002). 

Number counts and the redshift distribution are shown 
in Fig. 4. Note the sharp break in the QSO number counts 
at %'q ~ 19. The short-dashed line shows a power law with 
10 4m (which corresponds to a = 1 in Eq.[25]). As lens- 
ing magnification should be negative for shallower slopes, 
we cut the QSO sample to i'q < 18.8 (which contains 
952 QSOs for 0.8 < z < 2.5 within our mask) and use 
the fainter cut i'q > 18.8 data (which contains 671 QSOs 
within our mask) for comparison. Thus, unless stated oth- 
erwise, QSO are selected with 16.5 < i'q < 18.8, which for 
brevity will be sometimes refered as i'q < 18.8. 

Right panel in Fig. 4 shows how the redshift distribution 
for i' Q < 18.8 of EDR/QSO (histogram) and SDSS/QSO 
(dashed line) arc almost identical, with a mean redshift 
of ~zqso — 1-67. The redshift distribution for i'q < 19.1 
is very similar (with just higher surface density). Here 
we will only present results based on EDR/QSO. We have 
done all the analysis for SDSS/QSO and find identical re- 
sults in all cases. We choose to present EDR/QSO because 
it has been automatically produced (in the same pipeline 
as the galaxies) and because we are not concerned with the 
astrophysical nature of the QSO in the sample but rather 
with having a well defined photometric sample of distant 
objects. 

Figure 1 compares the galaxy pixel maps i' < 18.5 with 
the QSO i'q < 18.8 distribution in the same portion of 
the sky (QSOs arc the 2nd and 4th slices below the corre- 
sponding EDR/N and EDR/S galaxies). Figure 2 shows a 
zoom over the central region of EDR/N. Note that there 
is no apparent correlation between QSO and galaxies. As 
we will see below, most of the signal we are seeking for is 
hidden below the pixel resolution of this map. 

3. CORRECTING FOR GALACTIC EXTINCTION 

Absorption by Galactic dust can lead to a positive 
correlation between galaxies and high-z QSOs (see Nor- 
man & Williams 1999). Schlegel et al. (1998) extinction 
maps have a significant differential extinction E(B—V) ~ 
0.02 — 0.03 even at the poles. Thus, the extinction cor- 
rection has a large impact in the number counts for a fix 
magnitude range. The change can be roughly accounted 
for by shifting the mean magnitude ranges by the mean 
extinction. Despite this, extinction has little impact on 
clustering, at least for r' < 21 (see Scranton et al. 2001 
and also Tegmark et al. 1998). This is fortunate because of 
the uncertainties involved in making the extinction maps 
and its calibration. Moreover, the Schlegel et al. (1998) 
extinction map only has a 6'.1 FWHM, which is much 
larger than the individual galaxies and QSO we are inter- 
ested on. Many dusty regions have filamentary structure 
(with a fractal pattern) and large fluctuations in extinc- 
tion from point to point (ie see Fig. 3). One would expect 
similar fluctuations on smaller (galaxy size) scales, which 
introduces further uncertainties to individual corrections. 
As we are looking for a very low signal (of order of 1%) we 
have to be very careful with small systematic effects, such 
as extinction. 

The net effect of extinction is to produce a magnitude 

3 A similar argument and analysis can be made for other systcmati< 



absorption A which translates into a change in the local 
galaxy surface number density N: 

Ai^^A (1) 

where a is the slope of the number density counts: N ~ 
lQ am as a function of magnitude to. If the absorption 
A were a constant, this will would only change the mean 
number density in an uniform way. Unfortunately extinc- 
tion is highly variable and produces density fluctuations 
in the sky: 

S A = A A - < A A >=a(A- < A A >) (2) 

We thus see that extinction will introduce spurious num- 
ber density fluctuations in both the galaxy and the QSO 
distributions: 

~5~g = 5 G + s G 8 A (3) 
~5q=Sq + s q 5 a (4) 

where S G and Sq stand for the total observed fluctuations 
(as opposed to the intrinsic ones, Sq and Sq) and s G and 
sq are constant numbers characteristic of each population. 
From the above analysis we can see that extinction will 
introduce artificial cross-correlations between the galaxy 
and QSO populations, even when they are intrinsically 
uncorrelated. 3 One can correct for this type of effects by 
using the absorption maps. We can calculate the cross- 
correlations: 

<S G S A >= s g <S 2 a > (5) 

< S Q 5 a > = s Q <5\> (6) 

< 5 G S Q > = < 5 G 5 Q > +s G s Q <5\> (7) 

where we have assumed that the intrinsic galaxy and QSO 
positions are uncorrelated with extinction. We thus have: 

<SgSq>=<SgSq>+ < ~ Sg6a /<1 A ~ Sq> (8) 

where < 5 G Sq > is the intrinsic cross-correlation (eg from 
cosmic magnification). As we can measure all < 5 G 5a >, 

< Sq 5a > and < 5\ > from the maps, the above expres- 
sion allow us to correct the galaxy-QSO cross-correlation 
for extinction. 

To test this model we will consider 3 types of magnitude 
corections: raw magnitudes (ie uncorrected from extinc- 
tion) : i'\ extinction A corrected magnitudes: i* = i' — A; 
and extinction over-corrected magnitudes: i' + A. The 
later case will give us a good diagnostic to test our method 
to estimate the contamination in the galaxy-QSO cross- 
correlation. 

The results for these quantities are shown in the left 
panel of Fig. 5. Note how < S\ > (shown as closed squares) 
is quite flat. This is because of the lack of resolution on 
scales 9 < 10'. On larger scales, Schlegel et al. (1998) 
find that the angular spectrum of their absorption maps 
fit well P(k) oc fc~ 2 - 5 at all scale which corresponds to 

< 5 A >oc 8 +0 - 5 , so that fluctuations grow only slowly with 
scale. 

i, such as seeing variations. 
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Fig. 3. — Pixel maps of galactic absorption (Schlegel et al. 1998) over the EDR/N (top two slices) and EDR/S region. Wc also show below 
each extinction map the EDR/SDSS mask of pixels with less than 0.2 mag extinction in i' and less than 2" seeing. 
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Fig. 5 shows the prediction in Eq.[8] as a continuous line, 
which is close to 2% for the top case (in all cases we are 
using i < 19 pixel maps for both QSO and galaxies). As 
can be seen in the figure, when we over-correct the mag- 
nitudes for extinction (top panel) the correction is larger 
than when we use raw (uncorrected) magnitudes (shown 
in the middle panel), as expected. This illustrates how 
extinction could induce cross-correlations in the QSO and 
galaxy maps. When we apply the right (negative) extinc- 
tion correction to the magnitudes (bottom left panel) we 
find a negligible correlation within the error-bars. This 
result agrees well with that of using the raw magnitudes 
with and extinction mask (ie masking out the high extinc- 
tion pixels, see below). Either of this two last options is 
acceptable for our analysis as the residual cross-correlation 
is negligible given the error-bars. 

In the right panel of Fig. 5 we show the galaxy-QSO 
cross-correlation wgq for maps with extinction over- 
corrected magnitudes i* < 19 (closed squares) with wgq 
for maps with raw i' < 19 (open squares). The dashed 
line shows again wgq with extinction over-corrected mag- 
nitudes after subtracting the correction in Eq. [8] (ie contin- 
uous line in the top left panel of Fig. 5). The agreement is 
quite good indicating that there is only marginal contam- 
ination of extinction in wgq when we use raw magnitudes 
(in agreement with the middle left panel of Fig. 5). 

These results gives us some confident that we can con- 
trol this type of contamination in wgq (and also provides 
a further test to our cross-correlation codes) . 

3.1. Extinction and seeing mask 

Following Scranton et al (2002), pixels (in 6' resolution) 
with larger mean seeing or larger mean extinction than 
some threshold value are masked out from our analysis. 
The final mask is the product of this seeing and extinction 
mask (ie shown in Fig. 3) by the sample boundary (shown 
in Fig. 2). We have tried different thresholds following the 
analysis of Scranton et al (2001). As a compromise be- 



tween precision and area covered, unless stated otherwise, 
we use 0.2 maximum extinction and seeing better than 
1.8' for r' < 22 and i' < 22. For r' < 19 we relax the 
seeing cut to 2 arc-sec to include more galaxies. For these 
brighter galaxies 2 arc-sec provides very good photometry 
for clustering analysis (see Scranton et al. 2001). 

Bottom left panel in Fig. 5 shows the cross-correlation re- 
sults after applying the extinction mask over raw i' < 19 
magnitudes. The resulting contamination in wgq is neg- 
ligible. 

4. CROSS-CORRELATION MEASUREMENTS 

For our statistical analysis, we will use moments of 
counts in cells, eg the variance: 

wgg(0) =< 6 2 G (9) > (9) 

where Sq = tig/ < n > —1 arc number density fluctua- 
tions on cells of size 9 (larger than the pixel map resolu- 
tion) and < n > is the mean number of galaxies in the 
cell. The average < ... > is over angular positions in the 
sky. We follow closely Gaztanaga (1994, see also Szapudi 
et al. 1995) and use the same software and estimators here 
for the SDSS. This software have been tested in different 
ways and the results confirmed by independent studies (eg 
see Szapudi & Gaztanaga 1998). 

For the cross-correlation we use (see also MBM02): 

w GQ (9) =< S G (0)S Q (6) > (10) 

Note that this is different from the 2-point cross- 
correlation: 1^2(^12) =< Sg{0i)Sq(8 2 ) >, where #12 = 
\&2 — In our case, both cells are at the same location 
in the sky, and the scale dependence comes from changing 
the cell size. In fact, wgq is just an area average over W2, 
and its amplitude at scale 9 is typically 20 — 30% higher 
that u>2 at 612 — (see Fig.l in Gaztanaga 1994). 

Also note that shot-noise cancels out for cross-correlation 
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Fig. 5. — Left Panel: The density variance in the extinction < <5^ > (closed squares), compared with galaxy-extinction cross-correlation 

< $G $ A > (open circles) and the QSO-extinction cross-correlation < 5q 8 a > (closed circles). The continuous line shows the contamination 
in the galaxy-QSO cross-correlation as predicted by Eq. [8]. The top, middle and bottom panels correspond to extinction over-corrected 
magnitudes, raw magnitudes and results for raw magnitudes with extinction and seeing mask (which agrees well with the extinction corrected 
result). Right Panel: galaxy-QSO cross-correlation wqq as a function of cell radius 6>, for objects selected with extinction over-corrected 
magnitudes i* < 19 (closed squares and continuous line) and raw magnitudes i' < 19 (open circles). The dashed line shows the prediction in 
Eq. [8] using < 8 A >, < 8 a 8 a > and < 8q 8 a > cross-correlations shown in the top left panel. 



wgq , but not for w GG and wqq where the variance needs 
to be shot-noise corrected (eg Gaztahaga 1994). 

The reason for using the variance w G q rather than the 
2-point 102(^12) function is twofold. Firstly one expects 
the variance to provide better signal-to-noise ratios when 
the correlation signal is small, just because wgq is an in- 
tegrated quantity. Secondly, Menard et al. (2002) have 
argued that wqq provides a better estimator to improve 
the accuracy of the comparison with the theoretical pre- 
dictions in the weak lensing regime. We have in fact per- 
formed a calculation of the 1^2(^12) cross-correlation over 
the same data, using the simple counts-in-cell estimator 
for W2(0u) (eg Eq.[36] in Gaztahaga 1994). The results 
were dominated by noise (sampling errors) and we have 
therefore decided to just concentrate our efforts on wqq- 

We will also measure the galaxy-galaxy-QSO 3rd order 
moment: 

wggq{0) =< S G (9)5 G (9)S Q (9) > (11) 

and the galaxy-galaxy variance around QSOs: 

w GG ;q{0) =< S G (9)S G (9) > QSO (12) 

The extra or excess variance is dchncd as (MBM02): 

A(0) = w GG . Q (6) - w GG {9) (13) 

which turns out to be a measure of the 3-point function 
w GGQ (0) ~ A(0) (Fry & Peebles 1980, MBM02). 

Errors are obtained from a variation of the jackknife 
error scheme proposed by Scranton et al (2001, Eq. [10]). 
The sample is divided into N separate regions on the sky, 
each of equal area (in our case they are cuts in RA such 
that the number of pixel cells is equal in each region) . The 



analysis is performed N times, each leaving a different re- 
gion out. These are called (jackknife) subsamples, which 
we label L = 1 . . . N. The estimated statistical covariance 
for w G q for scales 0^ and 9j is then: 

Covar{9 t , 9j) = < Aw GQ (9 z ) Aw GQ (9 3 ) > (14) 

N-~\ N 
= ^ r Y, Aw GQ&)Aw G > Q (0 1 ) 

L=l 

Aw L G Q{9 l ) = w L GQ {9 i )-w^Q{9 i ) (15) 

where w G q(9i) is the cross-correlation measure in the L-th 
subsample (L = 1 . . . N) and WGQ{9i) is the mean value 
for the N subsamples. The case i = j gives the error 
variance. Note how if we increase the number of regions 
N, then jackknife subsamples are larger and each term in 
the sum is smaller. We have checked that the resulting 
covariance gives a stable answer for different N. 

These errors has been shown to be reliable when tested 
against simulations in Zehavi et al (2002) and have the 
great advantage of being model independent. In our imple- 
mentation we use N = 20 (but we have also tried N = 10 
and N — 40) independent regions of equal size. 

4.1. Variance in the galaxy-QSO correlation 

The top-left panel of Figure 6 compares the measured 
values of w G q(9) for the different magnitude bins cuts 
r' = 19,20,21,22 with unmasked pixels with the results 
for r' < 22 and i' < 22 with the extinction and seeing 
mask. All cases are roughly consistent with each other. 
Differences will be studied in section §5, the point here is 
that there is a significant signal in all bands (very similar 
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Fig. 6. — Top Left panel: Galaxy-QSO cross-correlation wqq (symbols with error-bars) for different galaxy samples as a function of cell 
radius (which is slightly shifted in each case to avoid overlaps). Closed and opened circles correspond to i' < 22 and r' < 22 with pixels 
masked (less than 0.2 mag extinction and more than 2" seeing). Closed triangles, open triangles, closed squares and open squares, show the 
unmasked results for galaxies with limiting fluxes of i' = 22, 21, 20, 19 respectively. Bottom Left PANEL: wqq for r' = 22 galaxy sample 
and different redshift bins in the QSO sample: 0.8 < z < 1.5 (closed squares) and 1.5 < z < 2.5 (open triangles) and 0.8 < z < 2.5 (closed 
triangles). Top Right PANEL: Circles and squares with error-bars show galaxy-QSO cross-correlation wqq for r' < 19 (z ~ 0.15) and for 
i' < 22 (z ~ 0.35). Closed triangles show wqq in i' < 22 when we randomize the galaxy counts. For comparison we also show wqq = 0.28~ 9 
as a dashed line. Bottom PANEL: Covariance Covar(9i,02) for 9\ = 0.5' (continuous line) and 9\ = 1.0' (dashed line) as a function of $2- 



results are find for cuts in g'). Bottom left panel com- 
pares wgq(0) when we separate the QSO in high and low 
redshifts. 

The top-right panel of Figure 6 shows in more detail 
the results for r' < 19 (z ~ 0.15, squares) and for i' < 22 
(z ~ 0.35, circles). The later have been shifted in angular 
scale (up by ~ 2.2) to match the r' < 19 depth. Closed 
triangles show the i' < 22 results after randomizing the 
angular position of the galaxy counts-in-cells. This pro- 
duces a distribution with identical variance, skewness and 
higher order moments, but which should not be correlated 
to the QSO sample. Thus, we expect wqq = in this 
case, exactly as found within the error-bars. This provides 
a test for our code and statistics against systematics such 
as boundary problems in the comparison. Also provides 
an indication that the errors are not underestimated. 

The bottom-right panel in Figure 6 shows the covari- 
ance matrix Eq.[14] for B\ — 0.5 (continuous line) and 
6*i = 1.0 (dashed line) from the jackknife estimator with 
N — 40 (values for N = 10 are comparable but the scatter 
is larger). As shown in the Figure the dominant contribu- 
tion comes from the diagonal terms. This is because we 
are using well separated bins (each cell is 4 times larger 
than the previous one). We neglect the off-diagonal er- 
rors in a first interpretation of the data, but we note that 
a stronger covariance arises from bins 62 < 0\ than for 
82 > 6*i, as shown in the Figure for the 9\ — 1.0' case. 

All cases on the right panel of Figure 6 correspond to 
maps with extinction 0.2 and 2" seeing mask for r' < 19 
and 1.8" seeing mask for i' < 22 (and i'g < 18.8). 

Assuming no covariance, the significance of the detec- 
tion (against zero correlation) for the first 5 points in top- 
right panel of Figure 6 is about 4-sigma in r' < 19 and 



8-sigma in i' < 22. 

4.2. Faint and Bright QSOs 

Bottom left panel in Figure 7 shows a comparison of the 
results in the EDR/N and EDR/S slices with the combined 
sample. The agreement is good within the errors, with 
slightly stronger signal in EDR/S. This provides a test 
for different versions in the QSO selection. Runs 752/756 
(which correspond to EDR/N) and 94/125 (which corre- 
spond to EDR/S) have different target selection version: 
v2.2a and v2.7 respectively. 

Top left panel in Figure 7 shows how the masked i' < 22 
galaxies cross-correlated with QSO subsamples cut at dif- 
ferent i'g < 18.8 bands. All results agree well within the 
errors with stronger signal for the brighter QSO, as ex- 
pected from the steeper slope of the brighter QSOs (see 
Fig.4). 

Right panel in Figure 7 compares wqq (for masked 
i' < 22 galaxies) for our nominal 16.5 < i'g < 18.8 with the 
faintest QSOs in EDR/QSO: 18.8 < i' Q < 19.1. Note how 
the fainter QSO sample (closed circles) show no significant 
cross-correlation. This result is in fact expected (see §5 be- 
low) if the cross-correlation is truly due to weak-lensing as 
a ~ 1 for this faint sample, see Fig.4. 

4.3. Comparison to galaxy-galaxy variance 

Fig. 8 summarizes the main observational results in this 
paper. We compare wgq(&) (open circles) and wgg{8) 
(closed triangles) for bright i'g < 18.8 QSOs with faint 
(r' < 19) and bright i' < 22 galaxies samples. The dot- 
ted line shows (for comparison) the power-law wgq = 
0.19~ - 7 . The i' < 22 case is also shown as closed circles 
in the top-left panel of Fig. 6. 
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Fig. 7. — Top Left panel: galaxy-QSO cross-correlation wqq for V < 22 galaxies with different QSO sub-samples: 16.5 < i'q < 18.3 
(closed squares), 18.3 < i'q < 18.8 (closed triangles) and the combined 16.5 < i'q < 18.8 (open circles). Bottom Left PANEL: Comparison of 
the mean wqq (open circles) with the values in the North (closed squares) and South (closed triangles) EDR subsamples (in all cases i' < 22 
galaxies and ig < 18.8). Right panel: Galaxy-QSO cross-correlation wqq (symbols with error-bars) for i' < 22 galaxies with different QSO 

sub-samples: 16.5 < i'q < 18.8 (open circles) and 18.8 < i'q < 19.1 (closed circles). 
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All data is well fitted to power-laws taking into account 
the errors (and neglecting the covariance). We find: 



__ rflog WGG _ 

— til 



Pgg 



a — d\ogw G Q 
PGQ = dlog(e) 



-0.72 ±0.02 z~0.15 
-0.60 ±0.03 z~0.35 



-0.83 ±0.17 z~0.15 
-0.96 ±0.13 z~0.35 



(16) 



(17) 



4.4. The Skewness 



Fig. 9 shows the pseudo-skewness (symbols with error- 
bars): 

m - (is) 

W GQ\ 6 > 

and the normalized excess skewness (continuous lines): 



A(9) _ w GG . Q (0) - w GG {0) 



GQ 



(0) 



w 



GQ 



(0) 



(19) 



where w gg .q is the galaxy variance around QSOs. It is 
expected that both quantities, S' 3 and S3 , should roughly 
agree on large scales. In fact, when shot-noise can be ne- 
glected S3 = S' 3 - 1 (see MBM02). Within the errors this 
relation is in good agreement with Fig. 9. 
A fit of a constant skewness in Fig. 9 gives: 



S'z = 



20.6 ±5.7 
13.6 ±4.3 



0.15 
0.35 



(20) 



In Fig. 9 we have used the masked mapped with 0.2 red- 
dening and 2" seeing. Results for other choices of the mask 
parameters are similar. 



5. COMPARISON WITH PREDICTIONS 

Consider the case of power-law correlations: £(r) = 
(Tsl\ or P{k) ~ k n , with 7 = n + 3. In current mod- 
els of structure formation n or 7 vary only smoothly with 
scale, so this should be a good approximation if we limit 
our study to a small range of scales. We will normalized 
the amplitude as: 



a 0.2 



0.2Mpc/hy 



R 



(21) 



where (Tq 2 refers to the non-linear amplitude of mass fluc- 
tuations on scales of 0.2Mpc/h, which correspond to ~ 1' 
in the sky. 

Galaxies might not be fair traces of the underlaying mass 
distribution. This is in fact one of the main motivations to 
use weak lensing as a tool to study large scale structure. 
In the Appendix we give a prescription to parameterize 
the possible bias that relates galaxies and matter distribu- 
tions. For a power-law correlation, this parameterization 
allows for a shift 60.2 in the amplitude and also a shift 76 
in the slope of the galaxy correlation with respect to that 
of the matter distribution. 



5.1. Projection and lens magnification 

To simplify notation, and without lost of generality, we 
will give all expression for a flat universe, Q m + = 1, 
where the comoving angular distance r(x) equals the radial 
comoving distance x ( see Bcrnardeau, Van Waerbeke & 
Mellier 1997, Moessner & Jain 1998, for the more general 
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case). Also, by default, assume Q m — 0.3, in agreement 
with current observations (see §6). 

The x distance is given in terms of the redshift z by: 



dx 



dz 



E(z) 



c 



n m + (i + z)zn m (22) 



which can be used to map x = x( z ) an d z — z(x)- For 
z ~ 0.15 we find a mean x — 430 Mpc/h, while z ~ 0.35 
we have x — 960 Mpc/h. 

The projected galaxy fluctuation is: 



S G (0) 



d X W G ( X ) 6 G (X,0) 



(23) 



where Xh is the distance to the horizon and refers to 
either the angular position in the sky or the radius of a 
circular cell in the sky. In the later case, which is the 
one we will study, 8 is smoothed over a cell or radius 9 
(a cone in the sky). The galaxy selection function Wg{x) 
corresponds to the probability of including a galaxy in the 
survey and is normalized to unity: 



rXH 
Jo 



d X W G (X) 



-r 

Jo 



dz «g(z) = 1 



(24) 



where n(z) is the normalized redshift distribution. 

On the other hand, fluctuations in the flux limited QSO 
induced by weak lensing magnification, 5^, can be ex- 
pressed in terms of the weak lensing convergence JC: 



5^6) = 2 (a - 1) K{9) 
where a is the slope of the QSO number counts: 



N{> i' Q ) ~ 10 



0.4 a i. 



(25) 



(26) 



In our sample we find a least square fit of a ~ 1.88 ± 0.03 
in the range 18.4 < i'q < 18.9 (shown as continuous line in 
Fig. 4). 4 The convergence is given by a projection over the 
radial matter fluctuation S(x,0), that acts as a lens. This 
projection is an integral over the lensing magnification ef- 
ficiency, S(x), a geometrical factor ~ x ^~ x that depends 
on the QSO XQ radial distribution and fl m . We express 
this projection as: 



5^6)= d x £(x)S( X ,0) (27) 



where 



£( X ) = 2(a-l) 



3H Ct r , 



<l + z) X d X '±—±W Q ( X ') 
* L J x X 

(28) 

and Wq(x) is given by the normalized probability to in- 
clude a QSO in our sample: 



dx 



W Q {x) = f 
Jo 



dz uq{z) = 1 



(29) 



As noted in Bernardeau, Van Waerbeke & Mcllicr 
(1997), the crucial difference between galaxy projection 



in Eq.[23] and lensing magnification in Eq.[27] is that £ 
is not normalized to unity. Thus, besides the projection 
effect, which can be modeled as an stochastic selection 
function, we have an overall re-scaling of fluctuation am- 
plitudes (this has been used in Gaztanaga & Bernardeau 
1998 to produce simple non-linear weak-lensing simula- 
tions). Figure 10 compares £{z)/E(z) with uq{z) and 
uq(z) for the QSO and galaxy samples in our analysis. 
For QSOs the redshift distribution is the one measured in 
the EDR/QSO 16.5 < iq < 18.8 sample, while for galaxies 
we show the predictions in Dodelson et al. (2001). 

5.2. Mass-mass and galaxy- galaxy correlations 

Using Eq. [21] and within the small angle approximation( 
eg see §7.2.1 in Bernardeau et al. 2002), the variance of 
the projected moss fluctuations can be expressed as: 



w(6) =< 8 2 (6) >= al 2 A 6 1 ' 



where 



A = A( 7 ) = B( 7 ) T n 



W G 



(30) 



(31) 



= (3- 7 ) (4- 7 ) (6-7) r( 7 /2-l/2)r(l/2) 
{1> 5^ 2 3 -r 9 T( 7 /2) 

T 1 is a geometrical factor of order unity (T 7 ~ 0.7 — 0.8) 
that comes from the area average over the 2-point func- 
tion: 



T = 
1 i — 



r(5- 7 ) 



2tt 



x dx I dxj) (l + x 2 — 2x cos < 



and W G is given by the radial galaxy selection: 



W G 



dX W 2 G (X) X 1 -''D 2 (z) 



(32) 



(33) 



where D(z) accounts for the redshift evolution of the cor- 
relation function: eg in the linear regime D(z) is the lin- 
ear growth factor, in the stable clustering regime D 2 {z) = 
{l + z)-\ 

For the galaxy-galaxy variance we assume the bias in 
Eq.[A4], which allows for a shift 60.2 in the amplitude and 
also a shift 7^ in the slope of the galaxy correlation with 
respect to that of the matter distribution. Thus, we have: 

w GG (6) =< 8 2 G (9) >= b 2 , 2 a 2 , 2 Aqq fl 1 "^ (34) 

where Aqg is the same as A in Eq.[31] for the new slope: 
Agg — A(-y + 27b). 

5.3. QSO- QSO and galaxy-QSO correlations 

We next want to estimate the galaxy-QSO cross- 
correlation. We express the angular QSO fluctuations as: 



(35) 



where 



Sq(9) stands for the intrinsic fluctuations while 
are fluctuations induced by magnification bias. 



4 At fainter mi > 19 magnitudes the counts show a very sharp break which could be related to possible incompleteness in our sample. Thus, 
can not estimate the slope reliably for the faint QSO sub-sample 18.8 < i'g < 19.1. 
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Fig. 10. — The dotted line shows the predicted galaxy redshift distribution uq(z) (divided by 3 to be on scale) for z = 0.35 (right 
panel), which corresponds to r' < 22 (and also i' < 22) and for ~z = 0.15 (left panel), which corresponds to r' < 19. The continuous 
line shows the normalized QSO redshift distribution uq(z) as measured in EDR/QSO. Short-dashed line is the lensing magnification effi- 
ciency £(z) I E(z) for the shown uq(z). The long-dashed line corresponds to yj ri(; (z) £{z) / E(z) which roughly represents the efficiency of 
galaxy-QSO cross-correlation. 



Neglecting the QSO (source) and matter (lens) cross- 
correlation (which is negligible given the large radial sepa- 
rations), the observed QSO variance has two contributions: 

w QQ (9) =< 5 Q (8) S Q (9) >= u£ Q (0) + w^(9) (36) 
where: 

^< W >= <2 A QQ 9^ (37) 

with: 

A QQ = B( 7 ) T 7 I s (38) 

£^ = f XH d x £ 2 (x) X^D^z). 
Jo 

The above expressions can be used, for a known magni- 
fication, to separate the intrinsic from the apparent QSO 
clustering. 

The QSO and galaxy populations are well separated in 
radial distances so that we can neglect the intrinsic cross- 
correlation < 8g{&) Sq(8) >. We then have: 

w GQ (9) =< S G (8) 5 Q (9) >~< S G (9) Spiff) > (39) 
so that using Eq. [27] and assuming Eq.[A4] we find: 



wgq(0) = fro. 2 <7o.2 A, 



GQ 



0l — 7 — lb 



with: 



A GQ =B{ 7 + ^ b ) T 1+lb W G E 

WgT= [ XH d X W G ( X )£(x) X^-^D^z) 
Jo 



(40) 



(41) 



5.4. Measure of bias and matter fluctuations 

The combination galaxy-galaxy and galaxy-QSO cross- 
correlation will allow us to break the intrinsic degeneracy 
between biasing and matter fluctuations, eg between &0.2 
and Co. 2 as measured by galaxy surveys alone. Here, we 
propose to measure the four parameters that characterize 
bias and matter fluctuations in a narrow range of scales 
(the power-law approximation). These parameters are cr§ 2 
and 7 for the variance in non-linear mass fluctuations (eg 
Eq.[21]) and 60.2, 76 for the bias function as described in 
the Appendix, eg Eq.[A4]. 

We can estimate these four parameters, 00.2, &0.2, 7 and 
7b, from the angular observations of galaxy-QSO correla- 
tions in the following way. We first take the measured 
logarithmic slopes of w GG and w G q in Eq.[16]-[17] from a 
fit to the data in Fig. 8. These slopes and the used to find 
the intrinsic matter slope 7 and bias scale dependence 71,: 



7 = 1 + /3 GG - 2[3 GQ 
lb = Pgq - Pgg 

We can then obtained 60.2 and 00.2 as: 

w GG (9) A G q 

Do. 2 - 77y\ 1 V 

w GQ {9) A gg 

A GG 



(42) 



7 0.2 — 



_ W GQ\ 



wgg{9) A 2 gq 



07-1 



(43) 



If the power-law model is a good approximation these val- 
ues should not be a strong function of scale. This will 
provide a consistency test for our power-law approxima- 
tion. 

The above reconstruction scheme has been tested 
against mock angular simulations in Gaztanaga (1995). 
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Fig. 11. — In both panels: open circles correspond to r' < 19 (z = 0.15) and closed squares to i' < 22 (z = 0.35). Left Panel: Values of 
b/r de-projected from the ratio wgg / W GQ assuming a constant, but stochastic, bias. Dashed lines are the best fit constant value of b/r for 
each galaxy sample. The dotted line gives the best fit linear relation for both samples. Right Panel: Values of fro. 2 and (To. 2 de-projected 
using Eq. [43], as a function of scale. 



This method can be generalized to extract the shape 
from wgg{@) in the quasi-scale invariance approximation 
by performing local inversions around R ~ 6\, where \ 
is the mean depth of the sample (see Gaztahaga 1995 for 
details). We have also verified the validity of such approx- 
imations over the weak lensing simulations presented in 
Gaztahaga & Bernardeau (1998). 

5.5. Scale dependence bias 

Left panel of Fig. 11 shows b/r defined in Eq.[Al3] and 
de-projected from the data at each point as 



b _ w GG A GQ 
r 



WGQ A G G 



(44) 



ie assuming that b/r is scale independent (ie 7^ = 0). Each 
point is shown at scale corresponding to the mean depth 
in the sample x@- The resulting values of b/r show a ten- 
dency to increase with scale, which means that 7^ = is 
not such a good approximation. But note that this trend 
is not very significant given the large errorbars: a constant 
value of b/r is only ruled out with a 3-sigma significance. 

Right panel of Fig. 11 shows the recovered values of 00.2 
and 60.2 as a function of scale given the prescription in 
Eq.[43]. The recovered values are quite flat, in good agree- 
ment with the assumption that the correlation functions 
can be approximated locally by power-laws (eg Eq. [34] ) . 

5.6. Results for (To. 2; ^0.2; 7 ond ^ 
From Eq.[16]-[17] and Eq.[42]: 



lb 



1.94 ±0.34 z 

2.33 ±0.26 z 

-0.11 ±0.17 1 

-0.36 ±0.13 1 



0.15 
0.35 

- 0.15 

- 0.35 



(45) 



(46) 



While the mean values from right panel of Fig. 11 are: 
, J 0.09 ±0.02 z~0.15 U7 , 

0Q.2 — <, OOQ 4-0 09 ^~n^K v 4 '-' 



CO. 2 



412 ±97 z~0.15 
364 ±86 z~0.35 



(48) 



These values assume that clustering is fixed in comov- 
ing coordinates. Values for other clustering evolution, ie 
e ^ 0, can be obtained using in Eq. [33]: 



£ 2 (z)~(l + z)-( 3+e > 



(49) 



If we require the underlaying (z — 0) (T0.2 to be the same 
in both samples we find: 



e~ 0.8 ±1.3 



(50) 



which roughly agrees with stable clustering, but has a large 
error. 

6. DISCUSSION 

As we will show in some detailed here, the values we 
have recovered for CT0.2, &0.2, 7 and 7b appear in contra- 
diction with currently popular expectations. We should 
therefore discuss first the two more important systemat- 
ics in our analysis: error-bar estimation (section §6.1) and 
extinction (§6.2). 

In §6.3 we give a brief account of what are the currently 
popular expectations for the above observables. In §6.4 
we will show how the results we find in this paper seem 
to contradict the expectations. A possible solution for this 
contradiction is hinted in §6.5 and §6.6 as we consider clus- 
tering evolution in shape and amplitude. In subsection 
6.7 we discuss how our interpretation changes with fl m 
and find that our proposed solution can also explained the 
low value of the measured skewness. Finally in §6.8, we 
compare our results with other estimates. 
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6.1. Unrealistic or correlated error-bars? 

We have sliced the data, as displaced in Fig.l, in N = 
10 — 40 RA bins (horizontal direction in the Figure). Be- 
cause of the mask, this results in subsamples with different 
shapes and holes, which increases the sample-to-samplc 
scatter in the jackknife error estimation Eq[14]. Thus we 
believe that this approach is conservative and slightly over- 
estimates all the jackknife errorbars presented in the above 
section. We have tested this idea with the APM (Maddox 
et al. 1990) pixel maps. 

We compare the jackknife error over a given slice out 
of the APM map against the scatter from different slices 
(the APM can fit over 20 EDR slices). The jackknife error 
is comparable to the zone to zone scatter on scales larger 
than a few arc-minutes (where the mean density of galax- 
ies in a cell is larger than unity), but is about a factor 
of 3 larger than the zone to zone errors on the smallest 
scales (dominated by shot-noise fluctuations). We do not 
attempt to correct for this here, but rather take the conser- 
vative approach of using the jackknife errors. This means 
that in principle the significance of the analysis presented 
below can be improved with a more sophisticated treat- 
ment of errors (eg see Colombi, Szapudi & Szalay 1998). 

There is also an overall shift in the mean amplitude of 
the zone to zone correlations (variance or skewness), due 
to large scale (sample size) density fluctuations, which is 
not taken into account by the jackknife error. The ampli- 
tude of this effect is comparable to the jackknife error on 
20' scales. This is similar to an overall calibration error as 
indeed corresponds to the uncertainty in the value of the 
mean density over the whole map (eg see Hui & Gaztanaga 
1999). Our error analysis hardly changes if we take this 
into account as we are mostly dominated by the (relatively 
large) scatter from shot-noise and small scale fluctuations. 

Other than this, the covariance shown in bottom-right 
panel of Figure 6 seems to be dominated by the diagonal 
terms and we have neglected it. 

6.2. Variable obscuration from the Galaxy 

Could our result be caused by small scale variations 
from obscuration in our galaxy? In §3 we have shown that 
large scale extinction is not affecting our results. There is 
a small but significant anti-correlation of the SDSS/EDR 
galaxy and QSO maps with Schlcgcl ct al. (1998) which in- 
duces an artificial galaxy-QSO cross-correlation. We have 
shown how this effect disappears when using an extinction 
mask. But Schlegel et al. (1998) maps have a FWHM of 
9 = 6', while we measure correlations below 8 — V . Could 
our cross-correlation be caused by extinction at smaller 
scales? The spectrum of extinction fluctuations as mea- 
sured by Schlegel et al. (1998), P(k) cx fc~ 2 ' 5 at scales 
from ~ 6 arc-min to ~ 30 degrees. , indicates that the 
amplitude of density fluctuations induced by absorption 
< 5\ >cx #+ - 5 grows slowly with scales. Thus extinction 
produces larger galaxy-QSO artifacts on larger scales. 

Extrapolating this spectrum to smaller scales should 
produce smaller contributions than the ones we have al- 
ready estimated in §3. 5 In summary, the detected galaxy- 
QSO cross-correlation have therefore the wrong amplitude 

5 In a recent paper, Kiss et al. (2002) find galactic cirrus emission 
stippcr spectrum at smaller scales. The extreme case a ~ 2 produces 
for our observations. 



and the wrong shape to be explained by galactic extinc- 
tion. 

6.3. The standard ACDM picture 

In the non-linear regime dark-matter (ACDM) numeri- 
cal simulations typically find that on cluster scales (r ~ 0.1 
Mpc/h) the non-linear variance of density fluctuations is 
— 10 2 £ L , where £ L stands for the linear value (com- 
pare continuous line with short-dashed line in left panel of 
Fig.12). Moreover £ NL ~ R^ 1 ( ~ i?~ 3//2 at the non-linear 
transition) which is typically steeper than the linear values 
£ L - R- 1 ' 2 ( - R- 1 at the non- linear transition). This 
behavior is reproduced by the non-linear fitting formulaes 
(eg Peacock & Dodds 1996) and have been extensively used 
in weak-lensing predictions (eg Bartelmann & Schneider 
2001, Benitez et al. 2001). For galaxies, predictions vary 
from model to model but one can typically find that on 
non-linear scales (1 to 10 Mpc/h) £ G ~ £ L for blue (Late- 
type) galaxies and up to £ G ~ 10£jv L for the red (Early- 
type) population. This gives a range b ~ ^/Cg/Ctvl — 0.1 
and b ~ 3. This is also in rough agreement with the esti- 
mations in galaxy surveys. We will call this the standard 
ACDM picture for non- linear galaxy and mass evolution 
(see Cooray & Sheth 2002 for an excellent review on all 
the above ideas). 

In terms of the variance most of these features can be 
summarized in the left panel of Figure 12, which is a ex- 
panded version of Fig. 11 in Gaztanaga (1995). Note how 
APM galaxies seem to be anti-bias on small scales with 
respect the ACDM model, something that seems to be 
compatible with different models of galaxy formation (eg 
Jenkins et al. 1998, Seljak 2000, Scoccimarro et al. 2001a, 
Sheth et al. 2001 and references therein). 

6.4. Weak lensing magnification 

If we take the detected galaxy-QSO cross-correlation sig- 
nal as being purely due to weak lensing magnification, we 
find important challenges for the above picture, even in 
the orders of magnitude involved. These are summarized 
in right panel of Fig.12, which shows the recovered mass 
variance as a function of scale. We find b ~ 0.1 (which 
seems to agree with the above blue galaxy model above) 
but note that £ NL ~ 10 4 £ L (ie compared upper closed fig- 
ures, £atli with dashed line, £ L ) which is much larger than 
expected in the standard ACDM picture above. 

We also find some evolutionary trends. First of all, from 
Eq.[45] we find £ NL ~ R~ 2 which is steeper (3 sigma sig- 
nificance) than in the standard ACDM picture (slope of the 
dashed line in right panel of Fig.12). At higher redshifts 
(z ~ 0.35) the slope seems slightly steeper £ NL <~ R~ 23 , 
as expected if we assume that on average there are lower 
mass halos at higher redshifts, and therefore steeper pro- 
files (see Navarro, Frenk & White 1996). Obviously, the 
slope of the galaxy £ G distribution is the same as in the 
standard A CDM picture (as these are direct observations 
that we also reproduce in our analysis). 

The measured values of the amplitude today cto.2 — 400 
is about 100 times larger that expected in the standard 

to varie from field to field as P{k) <x k a with —5.3 < a < —2.1 and even 
a constant value of < 5^ > as a function of scale, which can not account 
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Fig. 12. — Left panel: Comparison of the linear (continuous line) and non-linear (open squares) variance in the ACDM simulations 
(5 ) with that in APM galaxies (£ PM, squares with error-bars), from Fig. 11 in Gaztanaga (1995). Dashed lines shows the non-linear 

fitting formulae of Smith et al. (2002). The inset shows the bias function b(R) = y j^^ r for the non-linear (closed squares with errorbars) 

and linear ACDM (continuous line) cases. Right PANEL: Comparison of the linear (continuous line) and non-linear (dashed line) variance in 
the ACDM model with that from the bj < 20 APM galaxies (open squares with error-bars) and r' < 19 SDSS galaxies ( open circles with 
error-bars). Upper closed circles and triangles show the non-linear variance reconstructed from the galaxy-QSO correlation in the EDR/SDSS 
at z ~ 0.15 and z ~ 0.35 (scaled to z = under stable clustering). Lower closed figures show the corresponding results for the bias function 
b(R). The dotted line shows the non-linear mass variance prediction for a steep n e ff ~ 0.97 initial spectrum. 



A CDM picture, although the discrepancy factor depends 
on the particular cosmology. The significance of this result 
is over 4-sigma. 

The biasing amplitude b at 0.2 Mpc/h remains constant 
b ~ 0.1 with redshift. The null hypothesis of b = 1 is ruled 
out at a large significance. The slope -fb (ie from Eq.[46]) 
seems to become more negative at z ~ 0.35 (only 1-sigma 
effect), which partially masks the steepening of the density 
profiles as we increase the redshift. 

On weakly non-linear scales, the shape of the galaxy 3- 
point function and bispectrum can be used to determine b 
with independence of the underlaying spectrum. This was 
first proposed by Frieman & Gaztanaga (1994) and Fry 
(1994), following the ideas in Fry & Gaztanaga (1993). 
Recent measurements ( Frieman & Gaztanaga 1999, Scoc- 
cimarro et al. 2001b, Verde et. al 2002) find b ~ 1 at scales 
> 10 Mpc/h. In fact this value is not in contradiction with 
the values above as the bias seems to increase with scale, 
which extrapolates well to b ~ 1 at ~ 10 Mpc/h. 

Our recovered valued do not assume any model for the 
primordial spectrum or its subsequent evolution, but as- 
sumed a flat geometry with fl m = 0.3. This affects both 
the projection and the amplitude of the lensing magnifi- 
cation. The geometry enters in both the lensing and the 
galaxy projections and shouldjherefore cancel out to some 
extend for b, but not quite for ^ NL which is estimated from 
ratio Wqq/wgg- Thus we expect our estimations of b to 
scale with ~ -ffi and £ NL as ~ (fp-) 1 ' 5 ( see Bernardeau 
et al. 1997). Thus if fi m ~ 1 this would produce b ~ 0.3 
and (To. 2 — 158. These values are closer to the standard 
ACDM picture above but still uncomfortably far from it. 



Moreover they represent an inconsistency in the value of 

If the standard ACDM picture fails on these small scales, 
is there any possible alternative? 

6.5. Stable clustering 

The amplitude of £ NL at 0.2 Mpc/h follows the sta- 
ble clustering regime (structures are decoupled from the 
Hubble expansion): £ ~ (1 + z)~( 3+e \ with e ~ 0, indicat- 
ing that indeed clustering could be dominated by the halo 
profile (lh) rather than by the relative distribution of halo 
pairs (2h term in Eq.[86] of Cooray & Sheth 2002). The 
fact that the recovered amplitudes at redshifts ~z\ ~ 0.15 
and Z2 — 0.35 follow: 

M.fi±ny 3 .i. 6 (5i) 

does not involved any parameter fitting and is non trivial 
given that the raw observed correlation at different red- 
shifts are quite different (eg see Fig. [8]). In other words, 
there is no reason to expect any systematic effect to mimic 
this redshift dependence. Note that both in the linear 
regime £ £ ~ D 2 (; and in the strongly non- linear regime, 
this ratio should be independent of the amplitude £ n of 
the primordial spectrum. The fact that is good agreement 
with stable clustering on scales where clusters have col- 
lapsed is not surprising. Unfortunately the error we find 
for e in Eq.[50] is too large to draw more significant con- 
clusions. 

There is another aspect to stable clustering, which is 
related to the scale dependence of the correlations. Pee- 
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bles (1980) have shown that in the stable clustering for a 
power-law initial spectrum £ L <~ _R~(™+ 3 ) (where n is the 
index of the linear spectrum P{k) ~ k n ), the non-linear 
slope £ NL <~ i?~ 7 is: 



7 = 



3(n + 3) 
5 + n 



(52) 



In the ACDM models n ~ —1.5 on non- linear scales so 
that one would expect 7 ~ 1.3. Smith ct al. find that 
N-body simulations in fact fall a bit short of this predic- 
tion, with 7 ~ 1.0. For other values of n they find that 
7 increases with the steepness roughly as predicted by the 
stable clustering, but the measure 7 always falls a bit be- 
low the stable clustering prediction in Eq.[52] (sec Fig. 9 
in Smith et al. 2002). It should be noted that the agree- 
ment seems to improve for steep slopes: ie at n = 0, Smith 
et al. (2002) find 7 ~ 1.7 while Eq.[52] gives 7 = 1.8. Also 
note that Eq. [52] might be valid for stronger non-linearities 
(or smaller scales) than tested in the Smith et al (2002) 
simulations. 

Thus, the value we find from the galaxy-QSO data, 
7 ~ 2 in Eq.45], indicates a spectrum n ~ 1 in Eq.[52]. 
This value is very different from the standard A CDM pic- 
ture. But we will now show that in fact, for such steep 
initial spectrum normalized to as — 1, one expects the 
non-linear amplitude to be (T0.2 — 400, just as found in 
the our interpretation of the galaxy-QSO correlations. 

6.6. Non-linearities and halo model 

One of the possible origins of the discrepancy that we 
find for the measured value of (Jo. 2 as compared with the 
standard A CDM picture is the modeling of the non- linear 
variance £ NL . Based in the stable clustering ideas, Hamil- 
ton et al (1991) proposed that the non-linear variance 
should be a fixed universal function Jnl of the linear vari- 
ance once we identify the adequate mapping of linear to 
non-linear scales: £ NL ~ /jvl(£l)- Under stable clus- 
tering, the cosmic evolution of £ NL with the scale fac- 
tor a is £ NL <~ a 3 while £ L <~ a 2 , which implies that 

— —3/2 

£ NL ~ £ L . These transformations were generalized to 
the power spectrum and to models with ^ 1 by Peacock 
& Dodds (1994, 1996) and Jain, Mo & White (1995) to 
provide phenomenological fitting formulaes for non-linear 
clustering which have since been widely used. The range 
of validity and the accuracy of these initial ansatze was 
limited and, when looked in detail, there appears to be a 
fundamental problem in getting a universal fit that works 
for scale dependent models (see Baugh & Gaztahaga 1996, 
Smith et al. 2002). 

A more recent realization of this idea involves the dark 
matter halo approach (see Cooray & Sheth 2002 for a re- 
view) which assumes that the density field can be modeled 
by a distribution of clumps of matter (the halos) with some 
density profile. The large-scale clustering of the mass in 
the linear (and weakly non-linear) regime is given by that 
of the correlation between halos, which should trace the 
linear (and weakly non-linear) perturbation theory predic- 
tions, while the small-scale clustering of the mass arises 
from a convolution of the halo density profile with itself. 
Thus, in this model the results from weak-lensing mag- 



nification in right panel of Fig. 12 mostly traces the halo 
profiles. 

Smith et al. (2002) present new versions of the univer- 
sal non-linear fitting formulae Jnl that incorporate these 
new ideas. The new formulae seems to perform quite well 
against the Virgo N-body simulations. For the ACDM 
model, the comparison involves scales up to k ~ 20 h/Mpc, 
which corresponds to R ~ n/k ~ 0.16 Mpc/h, comparable 
to the ones of interest here. We will therefore use this new 
formulae for //vx (as implemented in the software made 
publicly available in Smith et al. 2002), to compare to our 
results. As a first test, left panel of Fig. 12 shows the Smith 
et al. (2002) non-linear prediction (dashcd-line) against 
our N-body results from ACDM simulation (opened cir- 
cles). The agreement is excellent and does not involve 
any additional parameter fitting. It should be noted that 
the differences with previous fitting formulaes (eg Peacok 
& Dodds 1996) in this particular case are small. This is 
important as it means that ACDM predictions for galaxy- 
QSO cross-correlations on non-linear scales in the liter- 
ature are correct and we need to look elsewhere for the 
origin of the discrepancy. 

As can be seen in the right panel of Fig. 12 the weak- 
lensing magnification predictions show a very steep non- 
linear variance, 7 ~ 2, which in the context of stable model 
or direct simulation (see §6.5 above) can be interpreted 
as indication of a steep linear spectrum n e // — 0.9 — 1.0. 
This steep linear spectrum can be obtained from the Smith 
et al. (2002) formulae by using a large effective value of 
the shape parameter T. For example, T = 200 (as = 1, 
TO = 0.3 and Oa = 0.7) produces 6 n e // = 0.96. This 
prediction is shown as a dotted-line in left panel of Fig. 12. 
This illustrates the point that we need a large positive 
slope to fit the data. Note that, once we fix the slope 
(from the same measurements) the only free parameter 
is ag, which we take to be as = 1 for consistency with 
power spectrum measurements (eg Gaztanaga & Baugh 
1998, Tcgmark & Zaldarriaga 2002 and references therein) 
and with the higher order correlations (see Bernardeau 
et al. 2002 and references therein). 

This means that the shape and amplitude for £ NL that 
we find from weak-lensing are compatible with the gen- 
eral picture of non-linear evolution of clustering and with 
the as = 1 normalization. Under this assumptions the 
only difference with the standard ACDM picture is that 
the linear spectrum on small (non-linear) scales has to be 
steeper. If we match this with the standard ACDM P(k) 
on larger scales, this requires P(k) to turn around and in- 
crease (rather than decrease) on galactic scales. As these 
scales are not strongly constraint by other observations 
(see below), the contradiction can be resolved in this way. 

6.7. Q m dependence and Skewness 

Other alternatives to the above considerations on non- 
linear modeling are to change the value of £l m and Q\ but 
keeping n e ff ~ — 1.5 (eg with a low shape parameter). 

Open models with Q m ~ 0.01 have more power on small 
scales, just as we need to explain the data. This does not 
seem to work for two reasons. First, because if we lower 
weak lensing magnification signal becomes smaller 
which translates into a higher value of the recovered mass 



"This large fiducial value for V illustrates how different is the measured slope from the expectations in the standard CDM picture. 
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£jv£. Second, because the non- linear slope 7 is too low 
compare to the weak-lensing reconstruction. 

As mentioned above (end of §6.4) closed models with 
Q m ~ 1 predict a three times stronger weak lensing magni- 
fication signal. This translates into a nine times smaller re- 
constructed variance and three times larger biasing, which 
is still not able to reconcile the discrepancy with the stan- 
dard predictions for the amplitude of non-linear matter 
fluctuations. 

Bernardeau et al (1997) predicted S3 ~ 
— 40S1~°- 8 zqso ~ 1,35 for the convergence skewness S3 in 
the ACDM model. One should roughly expect S' 3 ~ —S3, 
so that a direct comparison with Eq.[20] using our mean 

zqso = 1-67 yields Q m ~ 1.08 +o!65- Tnis is high 
and about 3 sigma out of the currently favored value 
of Cl m ~ 0.2 — 0.3 (eg Tegmark & Zaldarriaga 2002 and 
references therein). Note nevertheless that this is just a 
very rough comparison, both because we need to take into 
account non-linear effects (see Hui 1999) and because of 
the differences with pseudo-skewness. A direct comparison 
can be made to Fig. 3 in Menard, Bartelmann & Mcllicr 
(2002), which corresponds to zq ~ 1.5. Again here, the 
values for S3 in Eq. [20] favor the Q m ~ 1 and seem about 
3-sigma away from the fl m ~ 0.3 value. The S3 predic- 
tions are insensitive to the linear (scale dependence) bias 
and the amplitude and shape of the linear mass spectrum. 
But note that they are quite sensitive to linear spectrum 
n e f f of fluctuations on non-linear scales as they relay on 
the so-called saturation value in hyper-extended perturba- 
tion theory (Scoccimarro & Frieman 1999). If we change 
the value from n e ff ~ —1.5 to n e // — 0.9, as suggested 
by the non-linear slope 7 (see §6.5 above) then we get a 
value for S3 which is about 5 times smaller. We thus have 
S 3 ~ 8f2~ a8 ZQ5 _1 - 35 which now yields: 

— 0.15 (53) 

in better agreement with the ACDM cosmology. 

6.8. Comparison with other results 

As mentioned in the introduction, previous measure- 
ments of galaxy-QSO associations typically find evidence 
for weak lensing magnification. In fact, our results seem to 
agree well with the ones presented in Benitez et al. (2001 
and references therein) which also seem at odds with cur- 
rent expectation (see Bartelmann & Schneider 2001). 

Are these results in contradiction with other measure- 
ments of the mass or of the bias? On these scales we 
have measurements of the spectrum from Ly-alpha forest 
(Croft et al. 2002), but they are mostly sensitive to the lin- 
ear regime, as they correspond to higher rcdshifts. Weak 
shear lensing results (Hoekstra et al. 2002) seem to agree 
well with the standard ACDM picture above, and therefore 
seems at odds with the results presented here. In particu- 
lar, compare b/r ~ 1 in Fig. 6 in Hoekstra et al. 2002 with 
our b/r ~ 0.1 in Fig. 11, which are both obtained under the 
same definitions of biasing parameters. Note nevertheless 
that our analysis is probing slightly smaller scales, and 
more importantly note that we do not assume any specific 
model for the underlaying matter fluctuations. If the value 
of ^nl is indeed as large as indicated by our analysis one 
would also expect a stronger aperture mass variance, but 
given the number of parameters involved in the compar- 
ison of data with theory, some other interpretations may 



still be possible (eg Croft & Metzler 2000, Mackey, White 
& Kamionkowski 2002). More work is needed to under- 
stand this apparent discrepancy. The values of S3 found 
by Bernardeau, Van Waerbeke & Mcllier (2002) from the 
VIRMOS-Descart data should also be compared with the 
pseudo-skewness we find in Fig. 9. 

7. CONCLUSION 

Our results for a strong galaxy-QSO cross-correlation 
passed several test presented in §3 and §4 and Fig. 6-7. In 
particular note right panel in Figure 7 which shows how 
the signal disappears for faint QSOs, as expected if pro- 
duced by weak lensing magnification. Also note in top-left 
panel of Fig. 7 that the brighter the QSO sample the larger 
the amplitude of the positive detection. 

Strong galaxy-QSO or cluster-QSO lensing magnifica- 
tion can also contribute to the positive cross-correlation. 
Multiple QSOs images or just strong magnification by 
chance alignment will typically give positive correlations, 
in addition to the magnification bias effect. The key issue 
here is the probability of strong lensing, which for a fixed 
QSO redshift increases with the number density of QSOs. 
Thus, one would expect this effect to be stronger for the 
fainter QSOs, as the mean density is larger and the red- 
shift distribution is very similar. This is contrary to what 
we find in Fig. 7. Moreover, Benitez et al. 2001 argue that 
for strong lensing wqq ~ fsuiGG, where fs is the fraction 
of QSOs in a sample that are strongly lensed. In Fig. 8 
we find wqq — wqg on the smallest scale for z — 0.35. 
This requires fs — 1, which is unrealistic given the small 
size of the Eisntein radius of galaxies or clusters. Non- 
linear corrections to the weak-lensing approximation are 
typically small for the variance (see Menard et al. 2002). 
Thus, weak lensing magnification seems the most plausible 
explanation for the detected galaxy-QSO correlation (for 
other interpretations see Burbidge et al. 1990). 

Fig. 8 and Fig. 12 summarized the main results pre- 
sented in this paper, which seems at odds with the stan- 
dard ACDM model, both because the large amplitude and 
the steep slope of the recovered variance on non-linear 
scales. Note that this broad interpretation is independent 
of the details in our modeling to recover the variance. As 
can be seen in Fig. 8 we find u>gq — wgg on the small- 
est scale for z = 0.35. This alone, indicates that the bias 
factor b has to be quite small, b << 1, as weak lensing is 
typically less efficient that unity. In §6 we argue that a pos- 
sible explanation for this strong correlation is that the ef- 
fective linear spectral slope, P(k) ~ k neff , on small scales 
is n e ff ~ 1, which reproduces the measured slope 7 ~ 2 
in agreement with Eq.[52]. This is much steeper than in 
ACDM, where n e ff — —1.5 and 7 ~ 1.3. Within currently 
accepted cosmological parameters ( cr 8 ~ 1, Cl m ~ 0.3 and 
fl\ ~ 0.7) adopting n e // — 1 can explain at the same time 
the slope and amplitude of the recovered variance in Fig. 12 
and the low skewness in Fig. 9. In the ACDM cosmology 
this corresponds to a bump or step on the power spectrum 
P(k) on galatic scales. A physical interpretation for this 
feature can be found in the early universe, where a bump 
or step in the primordial sprectum is expected (see Barriga 
et al. 2001). If such a primordial feature is really present, 
it would have important implications for models of galaxy 
formation. Our weak lensing analysis requieres the bias 
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b to be a strong function of scale decreasing from b ~ 1 
at large scales to b = 0.1 at galactic scales, thus hiding 
the presence of such primordial bump. Note nevertheless, 
that if biasing is more complicated than modeled here (eg 
Eq.A4), our quantitative interpretation could be different 
in detail (see the Appendix for a discussion). 

This is a preliminary analysis from a small fraction of 
the SDSS early commissioning data. To explain the de- 
tected excess galaxy-QSO correlation of ~ 10% with cor- 
related photometric errors, such us flat-fielding or scatter 
light, we would need very significant photometric varia- 
tions: RMS of at least Ai' ~ 0.3 on arc-minute scales 7 , 
which is larger than the nominal calibration uncertainties 
of 0.03 magnitudes (see Stoughton et al 2002). It is im- 
probable that such large photometric errors are present 
in the EDR/SDSS as they would have already been de- 
tected as an excess galaxy-galaxy correlation when com- 
pared to other surveys (eg see Gaztahaga 2002a, 2002b, 
Connolly et al. 2002, Scranton et al. 2002). Neverthe- 
less, the full SDSS catalog should be able to do a much 
better analysis of systematics that the one presented here 
and confirm or refute our findings with a high significance. 
A larger piece of the SDSS would allow a comparisons to 
larger scales, which is missing in our analysis because of 



the narrow width of the EDR strips. On larger scales, 
the weak-lensing magnification signal can be more directly 
compared to other indications of mass clustering, which 
will be an essential test for cosmology and for the weak 
lensing magnification nature of these cross-correlations. 
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APPENDIX 
GALAXY BIASING 

We will consider two different mathematical approaches to parameterize the effects of biasing, ie. how galaxy fluctuations 
S G trace the underlaying mass distribution S. More elaborated physical models are based on semi-analytical galaxy 
formation (see eg Baugh, Cole & Frenk 1996, Seljak 2000, Scoccimarro et al. 2001, Sheth et al. 2001, Cooray & Shcth 
2002) but the results are quite model dependent. 

The first model consists in assuming that bias is linear, but non-local: 

S G (x) = J d 3 x' b(x - x') S(x') (Al) 

Which in Fourier space is just: 

S G (k) = b(k) S(k) (A2) 
This translates into a scale-dependence bias on the power-spectrum: 

P G (k) = b 2 (k) P(k) (A3) 

For a power-law correlation P(k) and a power-law b{k) we have that the galaxy variance is 



where 60.2 and 75 characterize the amplitude and scale dependence of the effective bias: 

/0 2\ 7b 

b(R) = b . 2 [—) . (A5) 

If 7fc = this reduces to the standard linear biasing model £ G (R) = b 2 i(R). 
In the second model we assume that bias is non-linear, but local: 

S G (r) = F[S(r)]. (A6) 

For small fluctuations S < 1: 

Sa-hS + ^l J 2 + | 5 3 + 0[<5 4 ] (A7) 
which gives rise to (see Fry & Gaztanaga 1993): 

Z G = bU + cf + df + ... (A8) 

where c = (S^ki + b^b\ + O.5&2) depends both on the non-linear biasing parameters and £3, the reduced skewness of the 
mass. Thus, the local model formally gives 

C G = m (A9) 
which reduces to the standard linear biasing model £ G {r) = b 2 £(r) for £ — > 0. On non-linear scales we can write: 

l a *X? (A10) 
alog£ 

which reproduces Eq.[A4] with cj> = 2^ b and K = bl 2 /al^ +lb \ 

Thus we have shown that two quite different hypothesis about biasing drives to a similar expression, ie Eq. [A4] , at least 
in the power-law limit. It is therefore plausible to assume that a more generic non-local and non- linear biasing could also 
be cast with such a parameterization. We therefore adopt Eq.[A4] on the assumption that a power-law should be a good 
approximation for the limited range of scales we consider here. We will be able to somehow test this assumption with the 
data. 



19 



Stochastic bias 

The stochasticity of the galaxy- mass relation can also contribute to the observed galaxy clustering (see Pen 1998, 
Scherrer & Weinberg 1998, Tegmark & Peebles 1998, Dekel & Lahav 1999, Matsubara 1999, Bernardeau et al. 2002, and 
references therein), so that the relation between 5 G and 5 is not deterministic but rather stochastic, 

6 G (r)=F[6(r)]+e s (r) (All) 

where the random field e$ denotes the scatter in the biasing relation at a given <5 due to the fact that 5(r) does not 
completely determine 5 G {r). Under the assumption that the scatter is local, in the sense that the correlation functions of 
e«5 (r) vanish sufficiently fast at large separations (i.e. faster than the correlations in the density field), the deterministic 
bias results hold for the two-point correlation function in the large-scale limit ( Scherrer & Weinberg 1998). 
This stochasticity is usually characterized by a parameter r defined as: 

< SS G > 2 < 5 G 5 G > 

r = — — with b A = — — 

b < 55 > < 55 > 

which implies: 



r 



with b 2 = (A12) 

< 55 G > 



b < 5 G 5 G > 



(A13) 



Note that r and b are mean quantities and both contain information on the scale independence and the stochastic e. This 
provides a different parameterization of biasing than the one given in the above subsection and has been used by Hoekstra 
et al. (2002) to characterize the comparison of the weak-lensing and galaxy correlations. To provide a comparison we 
also present in Fig. 11 results for b/r define in such a way, but note that the meaning of b here is quite different from the 
one in Eq.[A2] (see Dekel & Lahav 1999). Note that the cross-correlation coefficient |r| < 1, which introduces important 
constraints: eg b/r < 0.1 implies b < 0.1. 

The parameterization presented in the above subsection in terms of &0.2 and 7f,, is not necessarily in contradiction with 
the idea that the stochasticity could play an important role in biasing. We choose & .2 and jb f° r simplicity and to keep 
the number of parameters to a minimum. At each scale we can recover from our analysis the above stochastic modeling 
by simply replacing our constraints for 60.2 by constraints b/r and our constraints for 00.2 by constraints on rero.2- As 
|r| < 1 we have that CT0.2 > which shows that our prediction for 00.2, assuming r = 1, is a lower bound. Thus the 

qualitative conclusion in §7, that er .2 is too large as compare to ACDM models, is unchange by this modeling of bias. 
Further work would be needed to disentangle more complicated biasing schemes. 



